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Abstract 

From a variational action with non-minimal coupling with a scalar field and classical scalar and 
fermionic interaction, cosmological field equations can be obtained. Imposing a FLRW metric 
the equations lead directly to a cosmological model consisting of two interacting fluids, where 
the scalar field fluid is interpreted as dark energy and the fermionic field fluid is interpreted as 
dark matter. Several cases were studied analytically and numerically. An important feature of 
the non-minimal coupling is that it allows crossing the barrier from a quintessence to phantom 
behavior. The insensitivity of the solutions to one of the parameters of the model permits it to 
find an almost analytical solution for the cosmological constant type of universe. 
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1 1. Introduction 

2 In this paper we have considered fermionic and scalar fields as sources of dark matter and 

3 dark energy, respectively. Two types of interaction have been modeled, an interaction between 

4 the scalar field and fermionic field of the Yukawa type and a non-minimal interaction between the 

5 scalar and gravitational field. Recent observations suggest a equation of state of the type a> < — 1 
e [Q]], which makes possible to consider cosmological models admitting dynamical equations of 
7 state allowing to cross the barrier a> = — 1. A minimal coupling to the scalar field is not enough 
a to achieve a barrier cross to the phantom zone, however a non-minimal coupling allows this 

9 feature. This type of coupling has been studied in inflationary scenarios and in Grand Unified 

10 Theories (GUT) |H[1|1. 

11 We have studied the implications that these interactions have, specially non-minimal cou- 

12 pling, regarding the cosmic evolution of the model. Two Ansatz have been used and explored for 

13 the energy density transmission which are proportional to each energy density. We have found 

14 solutions to the field equations in a FRW flat universe. For these solutions the bilinear scalar 

15 S = if/ip shows a monotonic decrease in an expanding universe, this is, proportional with a -3 . 

16 This fact does not imply spinor energy density to have the same behavior. By redefining the 
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energy densities, a new equivalent and symmetric system of equations is found which manifestly 
has a redefined conserved energy momentum tensor. Inevitably, these redefinitions introduce 
two critical values to the scalar field, which are impossible barriers for the field to cross over, 
therefore, the cosmological solutions are valid between these two values f^. 

For the non-minimal parameter £ + and the energy density transfer function Q proportional 
to fa, we found a de Sitter university evolution coming from a decelerated phase to an accelerated 
phase. There is no crossing of the barrier a> = — 1, except for a large interaction parameter A 
which doesn't have a clear physical ground. For the redefined energy density transfer function Q 
proportional to p^ phantom solutions are obtained, and when the scale factor is near a c the Hubble 
parameter and energy densities diverge, finding a future singularity. To consider a non-minimal 
coupling, allows crossing the barrier a> - -1, where there is a transient from a dark matter 
dominated universe to a dark energy dominated universe, quintessence, and then a phantom 
evolution. 

We have a future singularity, in which the Hubble parameter, redefined energy densities and 
state equations diverge for a finite size, this future singularity correspond to a type III ||6t]. 

In the first section we present the Lagrangian model and field equations derived from it, and 
are particularized for the metric FRW. In the second section the cosmological field equations 
are found from the field equations derived in the first section. Two cases are studied: minimal 
and non-minimal coupling. In the third section the cosmological solutions are analyzed for two 
types of Ansatz interacting function. In the last section a discussion of the results are presented 
along with the curves derived from the cosmological field equations. A quasi analytical solution 
for these systems is presented for which a more detailed derivation is attached at the Appendix 
section. 

2. Formalism and field equations 

In this section, a brief description of the techniques used to include fermionic and scalars 
sources in mutual interaction with the Einstein theory of gravitation are presented J^SUl- Due 
to the fact that the gauge group of General Relativity does not admit a spinor representation the 
tetrad formalism is invoked. Following the general covariance principle, a connection between 
the tetrad and the metric tensor g MV is established through the relation 

fcn, = e'ejifc*, a,b = 0,1,2,3 (1) 

where e" denotes the tetrad or "vierbein" and r\ a \, is the Minkowski metric tensor. Here and after, 
Latin indices refer to the local inertial frame whereas Greek indices to the bundle space At. The 
main objective of this work is to describe the behavior of fermions iff : At — > R 4 and scalars 
cf> : At — > K with self-interacting potential density V(jj/, iff) : At — > K and U((f>) : At — > R 
in presence of a gravitational field. The dynamics between fermions and scalars fields will be 
represented on the Lagrangian by a Yukawa type interaction through the map / (0) : At — » M. 
The action for this system is, 



S(g, if/, iff, (f>) = 




where 

X = -(1 - ^ 2 )R + i (0T%0 - V^Py) - mifriff + 

-V&, iff) + l -d^<p - U(cf>) - if,f(M, (3) 
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remarking that spinors are treated here as classical commuting fields and f 6 R is a parameter 
for a non-minimal coupling between gravitation and scalar fields iTToL Bll. Lagrangian density 
(0 uses natural units, i. e. %tiG = c = K = 1, also m is the bare fermionic mass, tfr = tf/^y 
and R denotes the scalar of curvature. The Dirac matrices can be generalized to curved space 
through the definition = e^y" satisfying {F\r v } = 2g /JV '. The covariant derivative is then 
V^i/r = df,tfr - £1^ and V^i/r = + tfrQ.^,, where the spin connection O^, is given by 

^ = -^[r^-4(^)]yV, (4) 

with V . denoting the Christoffel symbols. 

The field equations are obtained by varying the total action (f2]i with respect to the tetrad, 
spinor field and scalar field, respectively. By defining a = 1 - %<p 2 , the following equations are 
obtained: 



Rfiv ~ \siiv R = a %v, (5) 
3/(0) da (<p) 

D ^ + — +tfr — * = ~W R ' (6) 



i^v^-im+nm = (7) 

where T MV = 7* + T° + Tj% and 

Tftv = -d^dycp + ^g^dpcpff'cf) - g^Uicj)) - g^Uaicj)) + d^aicp), (9) 

and -Cd = 5 ( | Ar / ' v /i | A _ V^t^r^t/f) - - V(^, 0) is the Dirac Lagrangian. 
In the following, we consider a FLRW flat universe described by the metric 

ds 2 = dt 2 - aitfidx 2 + dy 2 + dz 2 ], (10) 

where a(f) denotes the cosmic scale factor. According to the metric dipt , the tetrad components 
read 

and Dirac matrices become 

r° = r° , r = -±-y, (12) 

a(f) 

from which the spin connection components are obtained, yielding 

Qo = , n,- = l -d{t) Tm , (13) 
3 



71 where a dot for the time derivative have been introduced. 

72 We consider now that the fields are homogenous and isotropic. This is based on the observa- 

73 tional fact that on a cosmological scale higher than 3QQMpc the fields appear to be independent of 

74 the spatial coordinates in a post inflation evolution [1]. In order to study this model exhaustively 

75 we will consider the cases of a minimal coupling £ = and a non-minimal coupling £ + 0. 

76 Fermion field equations dT) and © become 

&+\h$ = -*Xm + /(0))yV~*y ^, (14) 
2 

= *(m + /(«)^ + ifV , (15) 

77 where H = H(t) — ^ is the Hubble parameter. From equations ([i4l and (fl~5T > the following 

78 relation is obtained: 

| (#) + 3//(#) = /(^rV - <Ar ^) ■ (16) 

79 Note that considering a self interacting potential of the form V (tfr, iff) = V (0Fi/r) turn null the 
so right side of equation dT6b ; in fact, with this assumption equation ( fTol l can be immediately inte- 
81 grated, yielding 



(17) 

82 where 5 = ^ has been defined, 
as 3. Cosmological Field Equations 

84 We will consider two cases, a minimal coupling £ = (a — 1) and a non-minimal coupling 

85 £ * 0. 

86 3.7. Minimal coupling ^ = 

87 In a minimal coupling a = 1 which simplifies considerably the field equations d5]-[8]l. It is 
as known that the Einstein field equations (f5J) fulfill 

- = -i[p^+P^ + 3(p^+p^)], (18) 
a o 

3// 2 = p^+p^, (19) 

89 known as the acceleration (TTFt and Friedman's equations dI9V On the other hand, the scalar field 

90 equation © can be written 

+ 3ff0 = -^-<^r . (20) 

91 allowing the following identifications 

P0 = l -4> 2 + U(<p\ (21) 
Pl// = <A(m + /(0))iA+ V0A,«A), (22) 

= ^ 2 -U(4>), (23) 
l_dV 15V 

* = ^ + m*- V{ *^' (24) 
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From d2"TI-l2"4ii regarding the fact that V = -3HS ^ it is straightforward to obtain 

p + 3H(p 4 + p$ = -Q, (25) 
/V + 3H(pj, + p^) = Q, (26) 

where Q = ty-grxfrij) is recognized directly from the field equations as the interaction function. The 
equation of state for each field is defined by a> x = j-, where x — <P or X = <A- 

This last remark allows to interpret fermionic fields as sources of dark matter and bosonic 
fields as sources of dark energy. We will see however that this same identification can be made 
in the case of a non-minimal coupling. 

3.2. Non-minimal coupling £, + 

In the case of non-minimal coupling, assuming ai^, = 0, i.e. a dust type solution, equations 
d5]- [HJ take the following form 

Pt+3Hpt(l + a>t) = ^_ -Q, (27) 

p^ + 3H P ^ = Q, (28) 
Q = pS, )a -^, (29) 

equations d27l i and (1281 imply a non conservation of the energy momentum tensor 7), v , however it 
will be shown below that there exists a redefined energy-momentum tensor that is conserved. 
Equation d29l ) is easily obtained by considering the Yukawa type of interaction / (0) = /3(f) and 
equation ( TTTb . The last two equations can be combined to yield an exact differential in terms of 
the scale factor a 

P "--A - o, 

o / 



da \ f3S o 

-a 3 + c. (30) 



f3S 



where the integration constant c will be shown to be zero. 

As it can be seen, there exist a singularity point a = a c when 1 - £<p 2 (a c ) = or equivalently 
when Pl p{a c ) = ^rp- with £ > in which the fields and the gradient of the fields diverge, 
therefore the model is valid only for |0| < <p (a c ) or \(f>\ > 0(a c ). For £ > a restriction for the 
field equations and a limited class of solutions is obtained [5]. 

Symmetry of the last system of equations can be restored by the following procedure. In the 
non-minimal coupling case, equations ( [TBI and ( TT9l can be generalized to 

d 1 

- = ~g[P<A +Aa + 3 (AA + Pi/,J], (31) 
3H 2 = /3 +A,, (32) 

through the re-definition of the densities (|2TT i. ( l22l i and pressures ( 1231 . ( l24t to 



a~ L \-(f + U(4>) + 6€H(W\, (33) 



Ac ; 


-- uHti 

















(34) 



^(# + 2 +2H<f>(f>) j, (35) 

i ay _ \ 

(36) 



Note how the non-minimal parameter £ is more relevant for the equations related to the 
bosonic field d33l ) and d33b . Note also that we recover the minimal coupling equations d2TI-l24T> 
when £ = as it is expected. 

On the other hand, we have the following equation obtained by adding (1331 and ( f35l > and 
replacing d30T > 

p (l + (Of) = a- l {H 2 (f V (1 - 2# - 2f#V) - ^'a 2 ^-H\ (37) 

119 which will be used in numerical simulation to find the initial conditions for the scalar field (p. 

120 The equations of state are 

P(/> - (Ofpf, 

Pif, - <*typ<j, 

121 which yields the following redefined conservation laws 

P<t> + ' iH \P<f> + P<t>) = - 1 _ 2 p> ~ Q, (38) 
fi+ + 3H(fr+P*) = Y^P* + & (39) 

122 where Q = a^ x ^ij^i<\> has been defined. In order to simplify calculations a dust type of equation 

123 for the fermionic field will be assumed, i. e. lo^ — in equation d39t . 

124 One interesting point to observe is that, through these redefinitions the energy-momentum 

125 tensor is now conserved WT^y = 0, with = a~ l T^ lv . At this point it is clear that we can 

126 interpret fermionic fields as sources of dark matter and bosonic fields as sources of dark energy 

127 as claimed above. 

128 We see that positive acceleration imposes 

P0 + /V + 3<<W < 0, (40) 

129 and from this equation we find < —4 (1 + r), in order to accelerated expansion to make sense, 

130 where we have defined r = ^ called the coincidence parameter, with ro ~ I the actual value. 



131 4. Cosmological Solutions 

132 On the premise that minimal coupling can be achieve on non-minimal coupling equations 

133 imposing £ = 0, we will focus on solving the general case, assuming an Ansatz for Q. Let us 

134 consider the following cases found in the literature lfIll[l2l[l3l[l4T.ll5ll 

Q = lAHfo, (41) 

Q = 3AHfa, (42) 
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135 
136 
137 
138 



where A is a positive parameter, which means dark energy represented by the scalar field is 
transferring into dark matter fermions, as current observations suggest yj]. 

Note that a dust type of model i.e. a>$ — or equivalently = imposed on the equation 
i implies 

l(-dV dV 



139 in accordance with equation (JT7J, let us suppose that the self-interaction fermionic potential is 

140 of the form V(if/,ip) = V(if)i//). This symmetry seems to be natural considering the need of 

141 the system to interact between matter and anti-matter in a way that is insensitive to a charge 
u2 conjugation operation. This allows to write equation (l43l in the following way 

V(Wr) = V(S) = bS, (44) 

143 where b is an integration constant. 

u4 On the other hand by replacing this result in equation (l34l i and use equations (l30l > and d44l it 

145 follows 

V(S) = -mS -pcS, (45) 

146 where we can identify b = -m and c — 0. This is justified because the force derived from this 

147 potential must have a range of the order of ~ -j-. 

148 Remark 1. In all the numerical solutions of the following subsections the parameters are fixed to 

149 tL»0 o = -0.98, A = 0.001 and £ = 4 (conformal coupling). This is according to last astronomical 

150 observations 

151 4.1. Case I: Q — 3AHp,p 

152 Using d39l and ( |4TT > an analytical solution for p^(a) is found 

p l/ ,=Pf a i «- 1 \ (46) 

153 with 

4> = 0oa 3/i , (47) 

154 where we have identified cpo = ^ and ciq — 1. The A dependence arises as a direct consequence 

155 of the interaction term. Replacing these results on (f38b we obtain an equation for that can be 

156 solved numerically. 

157 From equation ( 137] ) an initial condition for <\> in terms of observational and the model param- 

158 eters ^, A is found to be 



00 = + 



(1 

(48) 



) (3^ 2 - 12£l 2 + 2£0 + 3f A (l + 6*7^ ) + f (1 + ^ 



A 



where p c = 3H^ is the actual density critical value at t — to. We note also that in this case the 



condition 1 - t;<f> 2 = produces a c = {^4>q) 6 ' > it follows 

<t> («c) = 



(49) 



161 Case I presents <f>o — 2.22 < cf> c . 




(a) 




Figure 1 



As can be seen, Figure [Tal shows the redefined densities (dashed line) and p^ (solid line) 
behave in a way that in the distant future only dark energy survives, i. e. scalar field. Figure 
[Tb] shows a de Sitter like universe evolution in a distant future. A transition from decelerated 



expansion to accelerated expansion occurs at a 
an approximately constant value of - 0.2. 



0.6, which is shown in Figure |2a] to go along 





Figure 2 

Under the conditions of the simulation, in Figure [2b] a quintessence type of solution is ob- 
tained for at > — 1 which is according to experimental facts. However, using ojq < -1 a phantom 
type of behavior is obtained which has no physical foundation. Figure [3a] shows a decreasing r, 
which is expected for a compete dark energy dominance in the future. 

Figure[3b]indicates a decreasing rate of change from dark matter to dark energy, to the point 
in which dark matter practically vanishes in the future when the source has totally transferred 
its energy. Figure |4a] shows that <p grows rapidly for a «: 1 growing then to a lesser rate of 
change than initially observed, according to the behavior of the densities p^ and p^, Q and the 
coincidence parameter r. In Figure l4b] the potential V (a) is constructed, the self-interaction is 
almost null for a » 1 . 

4.2. Case II: Q = 3 AH p^ 

In this case the solutions will be obtained numerically due to the nonlinearity of the system of 
differential equations. Similarly as for the Case I from ( l37l l the initial condition for <p is obtained. 





Figure 3 





Figure 4 



Case II presents two branches for possible initial conditions for 0o 



185 
186 
187 



191 
192 



= 3.45 > 
= 2.43 < 



The first branch has a very similar behavior to Case I and is a quintessence universe. The 
second alternative is more keen to study due to the fact that it has a phantom type of behavior 
and for parameters of the model allows the crossing to quintessence. 

In Figure[5a]the redefined densities p^ and p^ show a good behavior up to the neighborhood 
of the critical point a c where both solutions diverge. The overall tendency is that dark energy 
density is always bigger than the matter density. Recall that both solutions are valid only in the 
range of the second branch. In Figure [5b] it is clearly seen that there is a divergence due to the 
factor a -1 in the solutions. 

Figure [6a] shows that this solution in the remote past the universe was decelerated and then 
at a * 0.7 begins to accelerate with a great increase of its growing rate in the neighborhood of 
a c . Adding a non-minimal coupling favors the phantom evolution allowing the crossing of the 
barrier a> — -1, which is clearly shown in Figure[6bl as claim the recent observations |Q]]. 

Figure [7a] exhibits similar behavior to those of the previous solutions due to the presence 
of a singularity at <f) c , however in the observational range there is a tendency to stabilization and 
constancy of the solution. Note also that the maximum is observed to be at the same turning point 
of the acceleration a » 0.7. Figure[7b]presents the same divergent behavior in the neighborhood 




(a) (b) 

Figure 6 



197 of <p c and has a tendency to remain constant before the singularity. 

198 Figure [8a] shows that under general lines the behavior of the scalar field is the same as that 

199 in Case I. The self interaction potential in Figure [8b] shows a dual behavior with respect of the 

200 expansion parameter in such a way that it is attractive for a < 0.7 and repulsive for a > 0.7. It is 

201 worth to mention that the transition happens just in the transition from decelerated to accelerated 

202 universe. 

203 5. Discussion 

204 As it can be seen in Case I, minimal coupling does not allow to pass over the phantom barrier 

205 at = - 1 for the initial conditions chosen. This seems to be the general behavior, reinforced by the 

206 exploration of the parameter space for large ranges. As far as numerical simulations permit, the 

207 behavior for this type of model in the future is quintessence. In Case II, non-minimal coupling 

208 allows to pass over the cj = — 1 barrier for certain initial conditions, although the model seems 

209 to be very insensitive of the value of the coupling parameter £ (see Figure [9]) . This type of 

210 model delivers a type III cosmological singularity produced mainly because of the a factor. Both 

211 solutions found are according to the observational measures. 

212 The incorporation of the scalar field as source of dark energy and the fermionic field as 

213 source of dark matter seems to model very adequately the dynamics as well as the interaction in 

214 the cosmological models found. 
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Figure 7 





Figure 8 



Having said that, in Case II, the £ parameter does not take such a prominent effect on the 
energy density solutions. These behavior is shown in Figure [9] where the dependence of the 
redefined densities to the cosmological parameter a and the model parameter £ are plotted. This 
£ parameter insensibility can be take under account to find an almost-analytical solution in the 
case of cosmological constant universe. In the appendix is shown how to find a Green function 
that leads to the following almost-analytical solution for the scalar field 



<M f ) 



r 

I 



df 



a(t') 



(50) 



where 



(51) 



222 The importance of this solution is that being valid in the barrier u)<p = —\ permits the com- 

223 parison of solutions on either side of the barrier, this let us propose the following remark. 

224 Remark 2. Let the solutions of a cosmological model derived from the Lagrangian (0 be valid 

225 in the interval [a, b] c R, then a scalar field solution fulfilling 



(p{t) > fa (t) V t e [a,b], 
11 



(52) 



Figure 9 



226 belongs to a phantom type of universe, whereas a scalar field solution fulfilling 

<j>{t) <<f> e (f) Vte[a,b], (53) 

227 belongs to a quintessence type of universe 

228 To complete the view, an almost-analytical solution is obtained for the Hubble parameter in 

229 terms of the scalar function 



2 J , 



<Pe (f) 

dt', (54) 



J (&]{t')-\)a'{t) 

230 see appendix for details. 

231 6. Appendix 

232 We will show here how to find a general Green's function for a non-minimal coupling at the 

233 boundary of cosmological constant. Using equations 03) . (l35b . and the equations of state we 

234 have: 

P (l + = (j> 2 - 2{L H (\cj> 2 )) , (55) 

235 where the operator Lg = 4* + H (t) has been defined. 

236 Note that on the left hand of equation d55l) . as is positive defined then for a cosmological 

237 constant solution (a)$ = — 1) follows 

l^=L H (yi). (56) 
12 



The formal Green's function solution for the previous equation is 

i2 



if 





G(t,s)(-^] ds. 



(57) 



239 
240 



where the kernel G:R(g)R^>R fulfills £ H G (t, s) = 5 (t - s) and £ H = J| + H (t) | . 

Let us define the following auxiliary function K . R(^) R — > R: K(t,s) = JjG (f, s) which 
fulfills the following differential equation 



—K(t,s)+H(t)K(t, s) = 5(t-s), 
ot 

by performing a Fourier Transform equation d58l ) can be written 

d 1 

— (f, x) + H (f) (f, jc) = — — exp (-ixf) . 
ot V2^ 

The equation d59l is the Bernoulli's differential equation, which has a known solution 



(58) 



(59) 



K(t',x) 



1 



/exp J H(t"')dt" 



exp(-/f"x) dt" 



V2^ 



exp 



jH(t")dt" 



by using the Hubble's parameter and an inverse Fourier Transform it follows 

( 



a (s) 



2n a (f) ' 



(60) 



(61) 



245 then the kernel is finally found to be 



Git, s) = 



a is) f df 



r df 

J a(t>Y 



By replacing (l62l on ( l57l i 

t r 00 2 



(62) 



(63) 



Lo 

note that the expression in brackets only depends on the parameter £ which allows to write 



J ^(f0 



(64) 
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248 which is equation d50i l. where 9^ - J a (s) (-57) ds has been defined. 



249 Assuming 6^ is analytic then it must be of first order in £ or higher to ensure convergence, as 

250 is readily seen on equation ( f55l l. A Taylor expansion is therefore of the form: 



% = ^a,e (65) 



251 which is the equation (fSTb . 

252 On the other hand, equation d37l > can be written 

C* + + 3 T + (i_4( fl )) ^ + ^ ^ ^ = °' (66) 

253 which by aim of equation (l30t takes the form 

6////'(l-^|(a)) + 3^0f(fl) = O, (67) 

254 which finally can be written 



2 j , 



J (^( f0 _l) fl '( f )3 

255 which is equation d54i >. 



df, (68) 



256 7. Acknowledgements 

This research was supported by DIR01 . 1 1 037.334/201 1 PUCV and Fondecyt 1 1 10076 (SL), 

258 also DI 10-0009 of Direction de Investigation y Desarrollo, Universidad de La Frontera (FP), 

259 and also by DI1 1-0071 (YV). 



260 References 

261 References 

262 [1] E. Komatsu, et al., Seven- Year Wilkinson Microwave Anisotropy Probe (WMAP) Observations: Cosmological 

263 Interpretation, Astrophys. J. Suppl. 192 (2011) 18. arXiv: 1001 . 4538 doi : 10 . 1088/0067-0049/192/2/18 

264 [2] I, L. Buchbinder, S. D. Odintsov, I. L. Shapiro, Effective action in quantum gravityBristol, UK: IOP (1992) 413 p. 

265 [3] E. Elizalde, S. D. Odintsov, Renormalization group improved effective potential for finite grand uni- 

266 fied theories in curved space-time, Phys. Lett. B333 (1994) 331-336. arXiv:hep-th/9403132 

267 doi : 10 . 1016/0370-2693 (94) 90151-1 

268 [4] T. Muta, S. D. Odintsov, Model dependence of the nonminimal scalar graviton effective coupling constant in curved 

269 space-time, Mod. Phys. Lett. A6 (1991) 3641-3646. doi : 10 . 1142/S0217732391004206 

270 [5] V. Faraoni, Inflation and quintessence with nonminimal coupling, Phys. Rev. D62 (2000) 023504. 

271 |arXiv : gr- qc/0002091 doi : 10 . 1103/Phys RevD ■ 62 ■ 23504 

272 [6] S. Nojiri, S. D. Odintsov, Inhomogeneous equation of state of the universe: Phantom era, future sin- 

273 gularity and crossing the phantom barrier, Phys. Rev. D72 (2005) 023003. arXiv:hep-th/0505215 

274 |doi : 10 . 1 1 03/PhysR evD ■ 72 . 023003] 



14 



275 [7] A. de la Macorra, Interacting dark energy: Decay into fermions, Astropart. Phys. 28 (2007) 196-204. 

276 arXiv : astro-ph/0702239 doi : 10 . 1016/j . astropartphys . 2007 . 05 . 005 

277 [8] S. Micheletti, E. Abdalla, B. Wang, A Field Theory Model for Dark Matter and Dark Energy in Interaction, Phys. 

278 Rev. D79 (2009) 123506. arXiv: 0902. 0318 doi : 10 . 1103/PhysRevD. 79 . 123506 

279 [9] M. O. Ribas, F. P. Devecchi, G. M. Kremer, Fermions as sources of accelerated regimes in cosmology, Phys. Rev. 

280 D72(2005) 123502. [arXiv: gr-q c/0511099 doi : 10 . 1103/PhysRe vD . 72 . 123502] 

281 [10] K. Nozari, S. D. Sadatian, Late-time acceleration and Phantom Divide Line Crossing with Non-minimal 

282 Coupling and Lorentz Invariance Violation, Eur. Phys. J. C58 (2008) 499-510. arXiv: 0809. 4744 

283 |doi:10.1 140/epjc/sl0052-00 8-0767-3| 

284 [11] W. Zimdahl, D. Pavon, Interacting quintessence, Phys. Lett. B521 (2001) 133-138. arXiv:astro-ph/0105479 

285 |doi : 10 ■ 1016/S0370-2693 (01) 01 174- 1 1 

286 [12] L. P. Chimento, A. S. Jakubi, D. Pavon, W. Zimdahl, Interacting quintessence solution to the coincidence problem, 

287 Phys. Rev. D67 (2003) 083513. |arXiv : astro-ph/0303145 ||doi : 10 . 1 1 03/PhysRevD . 67 . 083513] 

288 [13] S. del Campo, R. Herrera, D. Pavon, Interacting models may be key to solve the cosmic coincidence problem, 

289 JCAP0901 (2009) 020. TarXiv: 0812 .2210l|doi : 10. 1088/1475 -7516/2009/01/020 

290 [14] N. Cruz, S. Lepe, F. Pena, Dark energy interacting with dark matter and a third fluid: Possible EoS for this compo- 

291 nent, Phys.Lett. B699 (2011) 135-140. |doi : 10 . 1016/j .physletb . 2011 ■ 03 . 049 1 

292 [15] N. Cruz, S. Lepe, F. Pena, Dark energy interacting with two fluids, Phys.Lett. B663 (2008) 338-341. 

293 [arXiv: 0804. 3777 doi : 10 . 1016/j .physletb . 2008 ■ 04 . 035 1 



15 



